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C/^ . Abstract 

i> : 

■ We analyze the topological and fermionic vacuum structure of four-dimensional 
i QCD on the lattice by means of correlators of fermionic observables and topological 

^ ' densities. We show the existence of strong local correlations between the topological 

■ charge and monopole density on the one side and the quark condensate, charge and 
, chiral density on the other side. Visualization of individual gauge configurations 

Q> I demonstrates that instantons (antiinstantons) carry positive (negative) chirality, 

■ whereas the quark charge density fluctuates in sign within instantons. 
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^ . 1 Introduction and Theory 

During tlie last two decades several models liave been developed to describe 
the basic properties of QCD, quark confinement and chiral symmetry break- 
ing. The most popular are the dual superconductor, leading to confinement, 
and the instanton liquid model, which explains chiral symmetry breaking and 
solves the Ua{1) problem [1]. Both models rely on the existence of topologi- 
cal excitations, monopoles and instantons. Instantons have integer topological 
charge Q which is associated with the homotopy group n^^SU^Nc)) = Z, in 
four-dimensional Euclidian space. It is related to the zero eigenvalues of the 
massless Dirac operator Ip via the Atiyah-Singer index theorem [2]: 

Q = Index(^ ) = n+ - n_ , (1) 
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where n+ (n_) is the number of zero-modes of positive (negative) chirahty. 
Apart from this famous connection of topology and fermionic degrees of free- 
dom, here we attempt to systematically shed light on the relationship between 
the sea-quark distribution and topological density. We do this by studying cor- 
relators of topological densities with fermionic observables of the form 'tpTip 
with r = 1,74,75. Those quantities are usually referred to as the quark con- 
densate, quark charge density, and the chiral density. 

Recently, it was demonstrated that monopole currents which constitute a dif- 
ferent topological excitation of compact SU(3) gauge theory, related to the 
homotopy group n2{SU{Nc)/U{l)'^'=^^) = appear preferably in the re- 

gions of non- vanishing topological charge density [3,4]. It has been discussed 
that both instantons and monopoles are related to chiral symmetry breaking 
[5-7]. In [8] it was shown that zero-modes do not appear in configurations 
which are not monopole dominated. Here, by interpreting the correlators of 
monopole densities and fermionic operators, we confirm the local interrelation 
of monopole and fermionic degrees of freedom. 

For the implementation of the topological charge on the lattice there exists 
no unique discretization. In this work we restrict ourselves to the so-called 
field theoretic definitions which approximate the topological charge in the 
continuum, 

in the following ways [9]: 

- ^ E ^ (3) 

with 

= U,.{x)U,^{x) , (4) 

for the plaquette prescription and 

Ol"la = U {x, fJ,)U{x+ii, u)U{x+fi+u, p)U{x+il+u+p, a) 

X U^(x+i>+p+a, ijL)U\x+p+a, u)U\x+a, p)U\x, cr), (5) 

for the hypercube prescription. We mention here that the topological charges 
employed are locally gauge invariant, whereas the monopole currents are not. 
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The lattice and continuum versions of the theory represent different renormal- 
ized quantum field theories, which differ by finite, non-negligible renormal- 
ization factors. A simple procedure to get rid of renormalization constants, 
while preserving physical information contained in lattice configurations, is 
the cooling method. The cooling procedure systematically reduces quantum 
fluctuations, and suppresses differences between the different definitions of 
the topological charge. In our investigation we have employed the so-called 
"Cabibbo-Marinari method" which consists in a local minimization of the 
gluonic action in the SU(2) subgroups of SU(3) [10]. 

In order to investigate monopole currents one has to project SU{N) onto its 
abelian degrees of freedom, such that an abelian U{1)^^^ theory remains [11]. 
This aim can be achieved by various gauge fixing procedures. We employ the 
so-called maximum abelian gauge which is most favorable for our purposes. For 
the definition of the monopole currents 'm{x, fj,) we use the standard method 
[12]. To extract abelian parallel transporters u{x,fi) after imposing the maxi- 
mum abelian gauge one has to perform the decomposition 

U(x,ii) ^c(x,ii)u{x,iJ,) , (6) 

with {N = 3) 

u{x,iJ,) = diag [ui{x,ii),U2{x,iJ,),U3{x,iJ,)] , 
Ui{x, jj) = exp [i arg Uii{x, jj) - ^i0(a;, /x) , 
(t>{x,ij) ^^^rgUu{x,ij) e(-7r,7r] . (7) 



Since the maximum abelian subgroup U{1)^~^ is compact, there exist topolog- 
ical excitations. These are color magnetic monopoles which have integer- valued 
magnetic currents on the hnks of the dual lattice: 

T^i{x,lj) = ^ XI argwj(n), (8) 

where denotes a product of abelian links Ui{x,fi) around a plaquette □ 

and f{x + fiiii) is an elementary cube perpendicular to the n direction with 
origin x + p,. The magnetic currents form closed loops on the dual lattice as 
a consequence of monopole current conservation. Prom the monopole currents 
we define the local monopole density as 

tJ.,i 
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The local quark condensate i/jiIj^x) is a diagonal element of the inverse of the 
fermionic matrix of the QCD action. Other fermionic operators are obtained by 
inserting the Euclidian 74 and 75 matrices. We compute correlation functions 
between two observables Oi{x) and 02{y) 

g(y -x)^ {0,{x)OM) - {Oi){0,). (10) 

Since topological objects with opposite sign are equally distributed, we corre- 
late the quark condensate with the square of the topological charge density, 
and similarly for the other quantities. 



2 Results 

Our simulations were performed for full SU(3) QCD on an 8^ x 4 lattice 
with (anti) periodic boundary conditions for the (fermionic) gluonic sector. 
Applying a standard Metropolis algorithm, we checked that tunneling between 
sectors of different topological charges occurs at reasonable rates. Dynamical 
quarks in Kogut-Susskind discretization with Uf — 3 flavors of degenerate mass 
ma = 0.1 were taken into account using the pseudofermionic method with 40 
updates for the production runs. We performed runs in the confinement phase 
at inverse gluon coupling /3 = 5.2, corresponding to a lattice spacing a ~ 0.2 
fm. This leads to a total extension of the lattice of 1.6 fm and an energy unit 
of 1 GeV. Measurements were taken on 2000 configurations separated by 50 
sweeps. 

Figure 1 shows results for the correlation functions of Eq. (10) with Oi a local 
fermionic observable (except in (d)) and O2 the topological charge density for 
different numbers of cooling sweeps [13]. The correlations of the two densi- 
ties arc given in lattice units which leads to the small absolute values. The 
screening masses obtained from fits of the correlators to an exponential are 
collected in Table 1. The fits were taken over the interval indicated by the lines 
in Fig. 1. They clearly hint at a linear relation in the logarithmic scale. The 
statistical errors of the standard deviation in the plots are generally small, 
leading to relatively small errors in the fit parameters. The extracted quanti- 
ties have to be interpreted as effective masses and reflect the nonperturbative 
gluon exchange. Although cooling of quantum fields is necessary to extract 
topological structure, the exponential slopes of the correlators between the 
topological charge density and both the local chiral condensate (a) and the 
absolute value of the quark charge density (b) are nearly cooling-independent. 
We find that the correlations of the color charge density ipip^x) and |'?/'^-(/'(x)| 
with the topological charge density are very similar, both in the slopes and 
the absolute values. This becomes clear because the quark condensate can be 
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Fig. 1. Correlation functions of the topological charge density with the quark con- 
densate (a), quark charge (b) and chiral density (c), and its autocorrelation (d). For 
clarity of the plots error bars extending to negative values have been omitted. In 
the plots (c) and (d) for zero cooling steps some points cannot be resolved or have 
very large errors. 

interpreted as the absolute value of the quark density. However, cooling (or 
some other kind of smoothing) is inevitable to obtain nontrivial correlations 
between the chiral density, Oi = ip'j^ilj^x), and the topological charge density 
(c). This can be expected since both quantities are correlated via the anomaly. 
The topological charge of a gauge field is related to the chiral density of the 
associated fermion field by 

Q = q{x)d'^x = m 'ip'y^'ijj{x)d'^x. (11) 



We have checked that this relation also holds approximately for the corre- 
sponding lattice observables on individual configurations. The autocorrelation 
function of the density of the topological charge < g(0)g(r) > (d) should be 
compared to < ip'~f5ip{0)q{r) > (c). If the classical t'Hooft instanton with size 
p is considered, the topological charge density is 

g(x) ocp^(x2 + p2)-4 . (12) 
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Table 1 

Screening masses in GeV from fits to exponential decays of the various correlators for 
several cooling steps. The fit ranges are indicated by the lines in the corresponding 
figures. 

On the other hand the corresponding density of fermionic quantities [1] 

ipip{x) oc ipj5ip{x) oc p^(a;^ + p^)~^ (13) 

is broader. This behavior is reflected in (c) and (d) where the corresponding 
correlators are compared. The local relation q{x) = mil)^i^il){x) does not hold 
as one might naively expect, the fermionic distributions are of longer range 
than the instanton profiles. We performed check runs with 400 fermionic up- 
dates but did not find significant changes in the correlators. 

In Fig. 2 we show the analogous situation to Fig. 1 with O2 being the monopole 
density. Screening masses from corresponding fits to exponentials are found in 
Table 1 favoring again an exponential decay of the correlations. As in the case 
with the topological charge correlators, the monopole correlators are nearly 
cooling-independent for the local chiral condensate (a) and the modulus of the 
quark charge density (b). This expresses the strong correlations between the 
monopole and the topological charge densities themselves [3] . The correlations 
of the chiral density with the monopole density (c) deserve cooling to resolve a 
signal. We again compare the topological charge - monopole density correlators 
(d) to the chiral density results (c). The behavior of their widths is compatible 
with Eqs. (12) and (13) described above to discuss Fig. 1 (c) and (d). 

We now turn to a direct visualization of fermionic densities and topological 
quantities on individual gauge fields rather than performing gauge averages. 
We pursue this in the following to get insight into the local interplay of topol- 
ogy with the sea-quark distribution. In a series of papers we found that at 
clusters of topological charge density, which are identified with instantons. 
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Fig. 2. Correlation functions of the monopole density with the quark condensate 
(a), quark charge (b) and chiral density (c), and with the topological charge density 
(d). In the plots (c) and (d) some points cannot be resolved or error bars below zero 
are not printed. 



there are monopole trajectories looping around in almost all cases for SU(2) 
and SU(3) gauge theories, with improved action, in the presence of dynami- 
cal quarks, even across the deconfinement phase transition [3]. By analyzing 
dozens of gluon and quark field configurations we obtained the following re- 
sults. The topological charge is hidden in quantum fluctuations and becomes 
visible by cooling of the gauge fields. For cooling steps no structure can be 
seen in q{x), the fermionic observables and the monopole currents, which does 
not mean the absence of correlations between them. After a few cooling steps 
clusters of nonzero topological charge density and quark fields are resolved. 
For more cooling steps both topological charge and quark fields begin to die 
out and eventually vanish. 

In Fig. 3 a typical topologically nontrivial configuration, consisting of an in- 
stanton and an antiinstanton, from SU(3) theory with dynamical quarks on 
the 8^ X 4 lattice in the confinement phase is shown after 15 cooling steps for 
fixed time slices. We display the positive/negative topological charge density 
by white/black dots if the absolute value exceeds certain minimal fiuctua- 
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Fig. 3. Aspects of the fermionic field for a fixed time slice in a spatial volume 
of (1.6 fm)^ for a gauge field with an instanton-antiinstanton pair drawn by white 
and black dots, (a) Tjjip{x) > (light grey dots) within the instantons. (b) 'ip^ip{x) 
alternates in sign (light and dark grey dots) in a single instanton. (c) ^^^'il){x) is 
positive/negative in instanton/antiinstanton. 

tions. Monopole currents are defined in the maximum Abelian projection and 
only one type is shown by lines. The left 3D plot includes the local chiral 
condensate ijjijj{x), indicated by light grey dots whenever a certain threshold 
is exceeded. One clearly sees that both the instanton and antiinstanton are 
surrounded by a cloud of ^tpipi^x) > [14]. The middle 3D plot exhibits the 
situation for the quark charge density ■ip'^ip^x) indicated by light and dark grey 
dots depending on the sign of the net color charge excess. One observes that 
ip'^iplx) alternates in sign already in one instanton implying trivial correla- 
tions < ip'^%lj{x)q{ii) >= (not shown in Fig. 1). The right 3D plot displays 
the chiral density ipj^ipix) again indicated by light and dark grey dots. One 
nicely sees that the positive instanton is always surrounded by a lump with 
Tp'-^^ipi^x) > and vice versa. 



3 Discussion 

Combining the finding of Fig. 1 and 2 showing that the correlation functions 
between fermionic and topological quantities are not very sensitive to cooling 
together with the 3D images in Fig. 3, we conclude that instantons go hand 
in hand with clusters of i/^Til^^x) 7^ , F = 1, 74, 75, also in the uncooled QCD 
vacuum. In summary, our calculations of correlation functions between topo- 
logical densities and the fermionic observables exhibit exponential behavior. 
Results for the condensate and the modulus of the quark charge correlators 
are almost identical, as expected, since the quark condensate reflects the abso- 
lute value of the quark charge density. These correlation functions show little 
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cooling dependence. 

The correlations unambiguously demonstrate that not only the local chi- 
ral condensate iljip{x) but also the quark charge ip^ip{x) and chiral density 
■075'0(x) take non- vanishing values predominantly in the regions of instantons 
and monopole loops. Note that for the chiral density this behavior is expected 
due to the anomaly. 

Visualization exhibited that the distribution of sea-quarks is drastically en- 
hanced around centers of nontrivial topology (instantons, monopoles) in Eu- 
clidian space-time. Our investigations of full QCD with dynamical quarks rely 
on lattices of spatial volume of (1.6 fm)^ with a relatively large lattice con- 
stant of 0.2 fm. The correlation functions were also computed in the off-axes 
directions and show little anisotropy effects. This restauration of rotational 
invariance is usually taken as a sign of set-in of asymptotic scaling [15]. Nev- 
ertheless, it must be emphasized that our results represent the situation on a 
finite lattice with finite quark mass without the extrapolation to the thermo- 
dynamic and chiral limit. However, since all such correlators turned out rather 
independent of the gauge group and choice of the action etc., we expect that 
they are generic. We are currently studying fermionic observables along in- 
dividual monopole trajectories, to further shed light on the puzzle how far 
monopoles are involved in chiral symmetry breaking. 
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